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Using the simple ETAS branching model of seismicity, 
which assumes that each earthquake can trigger other earth- 
quakes, we quantify the role played by the cascade of trig- 
gered seismicity in controlling the rate of aftershock decay 
as well as the overall level of seismicity in the presence of a 
constant external seismicity source. We show that, in this 
model, the fraction of earthquakes in the population that 
are aftershocks is equal to the fraction of aftershocks that 
are indirectly triggered and is given by the average number 
of triggered events per earthquake. Previous observations 
that a significant fraction of earthquakes are triggered earth- 
quakes therefore imply that most aftershocks are indirectly 
triggered by the mainshock. 



1. Introduction 

There is a growing awareness and an intense research ac- 
tivity based on the fact that a significant fraction of earth- 
quakes are events triggered (in part) by preceding events. 
In addition, a significant part of triggered events may be 
indirectly triggered by a previous event through a cascade 
process. What is then the relative role of earthquake interac- 
tions and triggering compared with the underlying tectonic 
driving forces? Is there a way to distinguish triggered earth- 
quakes from untriggered ones or to estimate the proportion 
of directly or indirectly triggered earthquakes? Here, we 
use the Epidemic- Type Aftershock Sequence (ETAS) model 
to offer a quantification of earthquake interactions. This 
model is based on the two best established empirical laws of 
seismicity, the Gutenberg-Richter and the Omori law. The 
ETAS model has been used in many studies to describe or 
predict the spatio-temporal distribution of seismicity and 
reproduces many properties of real seismicity (see [Ogata, 
1999] and [Helmstetter and Sornette, 2002] for reviews). The 
ETAS model assumes that the seismicity results from the 
sum of an external constant loading and from earthquakes 
triggered by these sources in direct lineage or through a cas- 
cade of generations. From this definition (see below), it is 
clear that the ETAS model is not only a model of after- 
shock sequences, as the acronym ETAS would make one to 
believe, but describes the global seismicity including back- 
ground and interacting triggered seismicity. We use this 
model to quantify (a) the fraction of triggered events rela- 
tive to the sources and (b) the fraction of indirectly triggered 
events with respect to the total triggered seismicity. 

Question (a) has been previously visited in order to pro- 
vide unambiguous definitions of aftershocks and to decluster 
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seismic catalogs. Several alternative algorithms for the def- 
inition of aftershocks have been proposed [see Molchan and 
Dmitrieva, 1992 for a review]. Gardner and Knopoff [1974] 
and Knopoff [2000] used a windowing method and found that 
2/3 of the events in the catalog of Southern California are af- 
tershocks. Reasenberg [1985] analyzed the central California 
catalog and found that 48% of the events belong to a seis- 
mic cluster. Davis and Frohlich [1991] used the ISC catalog 
and found that 30% of earthquakes belong to a cluster, of 
which 76% are aftershocks and 24% are foreshocks. Kagan 
[1991] estimated the ratio of dependent events in various 
catalogs (California and worldwide) using an inversion by 
the maximum likelihood method of the ETAS model. The 
proportion of dependent earthquakes of the first generation 
that he estimated displays huge fluctuations from 0.1% for 
deep events to 90%, but is often close to 20%. 

With respect to question (b), it has long been suggested 
that aftershocks may produce their own aftershocks, com- 
monly known as secondary or indirect aftershocks. The ob- 
servation of large and sudden changes of the seismicity rate 
after a mainshock [e.g. Correig et al., 1997] and the ex- 
istence of strong spatio-temporal clustering of aftershocks 
shows that a significant proportion of aftershocks may be 
triggered indirectly by the mainshock, that is, they may 
be aftershocks of aftershocks triggered by the mainshock 
[Felzer et al., 2003]. For instance in Southern California, 
the M — 6.5 Big-Bear earthquake occurred a few hours fol- 
lowing the Landers M — 7.3 event and has clearly triggered 
its own aftershock sequence. While each aftershock induces 
a negligible stress change by comparison to the mainshock, 
all aftershocks when taken together can significantly alter 
the stress field induced by the mainshock, so that most af- 
tershocks at large times after the mainshock are triggered by 
previous aftershocks of the mainshock. Felzer et al. [2002] 
estimated the rate of indirect aftershocks, from a compari- 
son of the Landers aftershock sequence with numerical sim- 
ulations of the ETAS model. They found that about 85% 
of the aftershocks of the Landers event were indirect after- 
shocks. This implies that the 1999 M w = 7.1 Hector Mine 
earthquake was triggered, not by the 1992 Mw ~ 7.3 Lan- 
ders earthquake itself [Felzer et al., 2002], but more likely 
by some of its direct and indirect aftershocks. Felzer et 
al. [2003] further analyzed the temporal evolution of the 
proportion of secondary aftershocks. They found that, af- 
ter a few days or weeks following a mainshock depending 
on mainshock magnitude, most aftershocks are secondary 
aftershocks. We now recall the formulation of the ETAS 
model and its main results on the importance of triggered 
seismicity. 



2. The ETAS model of triggered seismicity 

The present parametric form of the ETAS model used 
in this paper was formulated by Ogata [1988]. We refer to 
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[Ogata, 1999; Helmstetter and Sornette, 2002] for reviews on 
the ETAS model and for a discussion of the model parame- 
ters. The ETAS model assumes that a given event of mag- 
nitude rrii > mo occurring at time ti triggers other events 
in the time interval between t and t + dt at the rate 

0m, (t - U) = p(mi)®(t - ti). (1) 

$(t) is the direct Omori law normalized to 1 
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where c is a regularizing time scale that ensures that the 
seismicity rate remains finite close to the mainshock. The 
average number of aftershocks triggered directly by an event 
of magnitude m is 



p(m) = fclO Q(m - mo) , 



(3) 



where m is a lower bound magnitude below which no 
daughter is triggered. The model is complemented by as- 
suming that each earthquake has a magnitude indepen- 
dently chosen according to the density distribution P(m). 
The magnitude distribution is usually taken equal to the 
Gutenberg-Richter law P(m) ~ i()- b (™-™a) witn eventu- 
ally a cut-off for large magnitudes. The model can also 
be extended to include the spatial distribution of seismicity 
[Ogata 1999]. The key parameter of the ETAS model (1) 
is the average number (or "branching ratio") n of directly 
triggered earthquakes per mother-event. This average is per- 
formed over time and over all possible mother magnitudes. 
The branching ratio has a finite value for 9 > equal to 



dt / P(m)p(m)$(t)dm. 



(4) 



The normal regime corresponds to the subcritical case n < 1 
for which the seismicity rate decays after a mainshock to a 
constant level (in the case of a steady-state source). Note 
that the realized number of aftershocks for a given earth- 
quake is not n but depends on its magnitude, according to 
the function p(m) given by (3). 

The total seismicity rate (or intensity) X(t) at time t is 
given by the sum of the "external" source s(t) and of the 
aftershocks triggered by all previous events 



X(t) =s(t)+ <l>™i(t-U). 



(5) 
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This external source s(t) acts as an external driving force 
ensuring that the seismicity does not vanish. 

Taking the ensemble average of (5) over many possible re- 
alizations of the seismicity, we obtain the following equation 
for the first moment or statistical average N(t) of \(t) [Sor- 
nette and Sornette, 1999; Helmstetter and Sornette, 2002] 



N(t) = s(t) +n $(t - T)N(r)dT. 



(6) 



The average seismicity rate is the solution of this self- 
consistent integral equation, which embodies the fact that 
each event may start a sequence of events, which can them- 
selves trigger secondary events, and so on. 

The global rate of aftershocks including indirect after- 
shocks triggered by a mainshock of magnitude M occurring 
at t = is given by p(M)K(t)/n, where the renormalized 



Omori law K(t) is obtained as a solution of (6) with the 
general source term s(t) replaced by the Dirac function S(t). 
The solution for K(t) is given in [Helmstetter and Sornette, 
2002] and is illustrated in Figure 1. The effect of the cas- 
cade of direct, secondary, and later-generation aftershocks 
is to renormalize the bare Omori law $(£) ~ l/t 1+e into 
K(t) ~ 1/t 1 "" at early times t « t* where t* « cjl - n\' 1/e . 
The characteristic time t* is infinite for n = 1 and be- 
comes very small for n <C 1. Figure 1 also shows the rates 
Ni(t) of aftershocks of generation i, for i = 1 to 20. Tak- 
ing an ensemble average, we predict Ni(t) = p(M)$(i), 
N 2 (t) = /J n$(t - T)p(M)<f>(r)dT, and more generally 



Ni(t) 



^{t-r)N l - 1 {T)dT, 



(7) 



such that the total seismicity rate is reconstructed as the 
sum N(t) = X/i^i Ni(t). Figure 1 illustrates clearly the 
role and importance of the successive generation of indirect 
aftershocks in the construction of the global observable seis- 
micity. 

In real data, it is impossible to distinguish unambigu- 
ously aftershocks from background seismicity, or direct af- 
tershocks from indirect aftershocks. The distinction is only 
probabilistic. Each event results in part from the external 
loading and in part from the effect of all previous earth- 
quakes. Knowing the parameters of the model, we can how- 
ever estimate the probability that each event results from 
the external source or is an aftershock of a previous earth- 
quake [Kagan, 1991]. In the sequel, we estimate the ratio of 
triggered seismicity over total seismicity in section 3 and the 
proportion of secondary aftershocks over total aftershocks in 
section 4, and we show that these two quantities are equal 
to the branching ratio n. 

3. Proportion of aftershocks 

Let us consider the situation in which s(t) corresponds 
to a constant Poisson source process with intensity p, rep- 
resenting the effect of the external loading. Then, the ob- 
served seismicity results both from this constant source rate 
and from the direct and indirect aftershocks triggered by 
this constant external loading. In the regime n < 1, the 
global seismicity is stationary, with large fluctuations fol- 
lowing large earthquakes due to the triggered aftershock se- 
quences. The rate of aftershocks ro triggered directly by 
the tectonic source p is on average n = pn because each 
single event triggers on average n events, when averaging 
over all magnitudes. The rate of second generation after- 
shocks, triggered by aftershocks of the tectonic source, is 
r2 = rwi = pn 2 . At the i th generation, the rate of after- 
shocks triggered indirectly by the tectonic source p is given 
by Ti = pn 1 . Summing over all generations, the global rate 
Raft of direct and indirect aftershocks of the constant ex- 
ternal source in the sub-critical regime n < 1 is given by 



Raft. = X! r ' =/i XX = i3 

i=l i=l 



(8) 



The global seismicity rate R is given by the sum of the ex- 
ternal loading p and of the rate of aftershocks R a ft.'- 



R — p + Raft. — p + 



pn p 



1 — n 1 — n 



(9) 



The result (9) shows that the effect of the cascade of after- 
shocks of aftershocks and so on is to renormalize the exter- 
nal constant source p to a higher level R that increases as 
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n is close to the critical value 1, as illustrated in Figure (2). 
This result is well-known in the branching process literature 
[Harris, 1963] and has also been derived by Kagan [1991] for 
the slightly modified version of the ETAS model using c = 
and replacing it by an abrupt cut-off at early times. 

The proportion of aftershocks (of any generation) is thus 
equal to R a ft./R = n. This expression shows that the av- 
erage branching ratio n can be directly observed from a 
suitable analysis of seismicity catalogs. Indeed, clustering 
algorithms for detecting and counting aftershocks provide 
a direct estimation and in general a lower bound of n be- 
cause most triggered events cannot be distinguished from 
the background seismicity. Note that the result (9) can also 
be derived directly from the master equation (6) by inserting 
s(t) — (i in (6) and taking the expectation of N(t). 

4. Proportion of indirect aftershocks 

There is another interpretation for n as well as an addi- 
tional empirical tool to estimate it. We calculate the total 
number of aftershocks nt triggered by a mainshock of magni- 
tude M, including all the generations of direct and indirect 
aftershocks, as follows. The number of direct aftershocks 
is given by ni = p(M) using the definition (1). The aver- 
age number of second generation aftershocks 712 is given by 
the product of m with the average number of aftershocks 
per earthquake defined by n. Therefore ni = p(M)n. The 
number of third generation aftershocks of the mainshock is 
TI3 = p(M)n 2 . The number of aftershocks for the i gen- 
eration is m = p(M)n t ~ 1 . The total number of aftershocks 
triggered by a mainshock of magnitude M is thus given by 



S = Y,n i = P {M)Y,n i = P - 



p(M) 



(10) 



For k « 1, S R! p{M), i.e., most aftershocks are directly 
triggered by the mainshock. For n Ri 1, S > p{M), i.e., 
most aftershocks are indirect aftershocks of the mainshock. 
The proportion of indirect aftershocks is given by 



S-m 



KM) 
1-n 



p(M) 



P{M) 



(11) 



This result (11) shows the fraction among all aftershocks 
of the aftershocks triggered indirectly by the mainshock is 
given by the average branching ratio n, independently of the 
mainshock magnitude M. We can also derive the result (11) 
from the master equation (6). Inserting s(t) — S(t)p(M) in 
(6) and taking the integral of (6) gives after some manipu- 
lation the global number of direct and indirect aftershocks 
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N(t)dt = p(M) + n / N(t)<1t = p(M) + nS, 



which recovers expression (10) for S. 

The branching ratio n gives the proportion of indirect af- 
tershocks averaged over the whole aftershock sequence. It 
is different from the instantaneous proportion of indirect af- 
tershocks u(t) that is defined by 



"(*) = 



K(t) - $(t) 
K(t) ' 



(12) 



which can be computed analytically using the expression of 
K(t) given by Helmstetter and Sornette [2002]. The instan- 
taneous proportion of indirect aftershocks increases from 
for very small times (<c (all aftershocks are triggered di- 
rectly by the mainshock) to a maximum value smaller than 



one at large times t > t* given by 

u 00 =l^(t) = l-(l-n) 2 ^^. (13) 

The temporal evolution of v{t) given by (12) is illustrated 
in the inset of Figure 1. 



5. Conclusion 

We have shown that, in the ETAS model, the proportion 
of earthquakes that are triggered is equal to the proportion 
of aftershocks that are indirect, and is given by the branch- 
ing ratio. Previous observations that a significant fraction of 
earthquakes are triggered earthquakes therefore imply that 
most aftershocks are indirectly triggered by the mainshock. 
The importance of indirect aftershocks casts doubts on the 
relevance of prediction of aftershocks rate based on the cal- 
culation of the Coulomb stress change induced by the main- 
shock only, neglecting the stress changes induced by after- 
shocks [Stein, 1999]. It also opens the road for improved 
methods of seismicity forecasts [Felzer et al., 2003]. 
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Figure 1. A realization of the ETAS model showing the 
realized seismicity rate X(t) (circles) following a M = 7 
mainshock obtained by averaging over 1000 simulations 
with n = 0.8, a = 0.8, b = 1, 9 = 0.2, m = 0, c = 0.001 
day, and the average renormalized propagator K (£) (solid 
gray line). The bell-shaped curves show the seismicity 
rates Ni(t) of aftershocks of generation i estimated from 
equation (7), for i = 1 to 20 from top to bottom. The in- 
set gives the proportion of indirect aftershocks v(t) eval- 
uated by (12). After 15 minutes, most aftershocks are 
triggered indirectly by the mainshock. At large times 
Of, the proportion of indirect aftershocks goes to an 
asymptotic value of 0.97. 
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Figure 2. Rate of seismic activity for a synthetic cat- 
alog generated using the ETAS model with parameters 
/i = 0.1 source events per day, c = 0.001 day, n — 0.8, 
6 — 0.2, b = 1 and a — 0.8. The average seismicity rate 
is close to the expected value fi* = y,/(n— 1) predicted by 
(9) (dotted line) and is always significantly larger than 
the constant external rate fj, (dashed line) . 78% of earth- 
quakes are aftershocks, among which 79% are indirect 
aftershocks, in good agreement with the predictions (9) 
and (11). 
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